Mon. Not. R. Astron. Soc. 000, 000-000 (1993) 



A deep i^05'^r survey - IV. The evolution of X-ray-selected 
QSOs 

B.J.Boyle/ T. Shanks,^ I.Georgantopoulos,^ G.C. Stewart^ and R.E. Griffiths^ 

Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CBS OHA 
Department of Physics, University of Durham, South Road, Durham, DHl 3LE 
Department of Physics and Astronomy, The University, University Road, Leicester, LEI 7RH 
^ Department of Physics and Astronomy, Johns Hopkins University, Homewood Campus, Baltimore, MD 21218, USA 



Accepted 1994 July 7. Received 1994 June 27; in original form 1994 February 22 



ABSTRACT 

We report on a new estimate of the QSO X-ray luminosity function and its evo- 
lution with redshift based on a sample of 107 QSOs detected at faint X-ray fluxes, 
5(0.5 - 2keV) > 4x lO'i^ erg s'^cm'^, with the _R 05^ T X-ray satellite. Forgo = 0.5, 
the X-ray evolution of QSOs in this sample is consistent with strong luminosity evo- 
lution, L'^{z) x {1 + z)3 25±o.i^ at low redshifts (z < 1.60) and a constant comoving 
space density at higher redshifts. The derived rate of evolution at low redshifts is 
thus significantly higher than that obtained previously for the Emsiem Extended 
Medium Sensitivity Survey (EMSS). Indeed, most luminosity evolution models pro- 
vide a very poor fit (rejected at the 99 per cent confidence level) when applied to 
the combined EMSS and ROSAT samples, although a polynomial evolution model, 
_L^(z) cx )^ provides an adequate fit for go = 0. For go = 0.5, a simple 

power-law luminosity evolution model with a redshift cut-off {L%i^{z) cx (1 -|- 2;)2 5i±o 
Zjnax = 1.25) is an acceptable fit to the EMSS and _R 05^ T samples only if a sizeable 
dispersion in the QSO X-ray spectral index ((T(ax) ~ 0.5) is included. More complex 
evolutionary forms, such as a variable luminosity function slope or density evolu- 
tion combined with luminosity evolution, also fail to provide an adequate fit to the 
combined data-set. Possible systematic differences between the Emsiem and ROSAT 
energy bands (e.g. spectral index variations) may account for some of the observed 
discrepancy between the QSO samples. Based on the observed range in the parameter 
values for the best-fit evolutionary models, we obtain formal values of 34-53 per cent 
for the QSO contribution to the 1-2 keV X-ray background. 

Key words: X-rays: general - cosmology: diffuse radiation - quasars: general - 
galaxies: active 



1 INTRODUCTION 

The QSO X-ray luminosity function (XLF) and its evolution 
with redshift provide fundamental information on the X-ray 
properties of QSOs. There are now a significant number of 
deep X-ray pointings which have been obtained with ROSAT 
(Shanks et al. 1991; Hasinger et al. 1993; Branduardi- 
Raymont et al. 1993), and in this paper we present the 
results of analysis of the QSO XLF and its evolution based 
on a sample of 107 QSOs detected to an X-ray flux limit 
5(0.5 - 2keV) > 4 x 10"^® erg s"^cm"^ with the Position 
Sensitive Proportional Counter (PSPC) on board ROSAT 
(see Georgantopoulos et al. 1994). In an earlier paper in 
this series (Boyle et al. 1993, hereafter Paper I) we reported 
on an initial determination of the QSO XLF, $x(iX) ^)) 
based on a preliminary sample of 42 QSOs detected on two 



deep (~ 27 000 s) ROSAT fields with an X-ray flux limit 
of 5(0.5 - 2keV) > 8 x 10"^^ erg s'^cm"^. We confirmed 
that the XLF exhibits a two-power-law form with a 'break' 
luminosity at ix(0-3 - 3.5 keV) = 10'*^^^='=°^^ ergs"\ first 
identified by Maccacaro et al. (1991, hereafter M91) and 
Delia Ceca et al. (1992) from an analysis of 448 QSOs in 
the Einstein Extended Medium Sensitivity Survey (EMSS). 
We also confirmed that the evolution of the XLF was well 
represented by a power-law increase in luminosity with red- 
shift ix oc (1 -|- z)^ , but derived a higher rate for the com- 
bined EMSS//?05Ar samples (k = 2.8±0.1) than had been 
obtained by M91 (k = 2.56) based solely on an analysis of 
the EMSS. In addition, we found tentative evidence for a 
'cut-off' in the luminosity evolution of QSOs at z'i.2, but 
had insufficient QSOs at these redshifts to establish firmly 
its existence or nature. Uncertainties in the evolution of 
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QSOs at high redshift, coupled with the uncertainty in the 
value of the faint-end slope of the XLF, led to rather weak 
constraints (30-90 per cent) on the overall contribution of 
QSOs to the soft (2 keV) X-ray background. 

In an attempt to improve our knowledge of the XLF and 
its evolution (particularly at low luminosities and at high 
redshifts) we have recently increased the size of our sample 
to 107 QSOs (Georgantopoulos et al. 1994) by obtaining 
fibre-optic spectroscopic identification for the optical coun- 
terparts to faint X-ray sources detected in a further three 
ROSAT ^ieXAs, the deepest of which extends our X-ray flux 
limit to 5(0.5 -2keV) > 4 x 10"^® ergs"^cm"^. In this 
paper we present the results of a new analysis of the XLF 
based on this larger sample of X-ray-selected QSOs. In Sec- 
tion 2 we present brief details of the QSO sample used in 
this paper, and we describe our analysis and results in Sec- 
tion 3. We discuss these results in Section 4 and present our 
conclusions in Section 5. 



2 DATA 



2.1 ROSAT observations 

We base our analysis of the QSO XLF on a new sample 
of 107 QSOs detected with the PSPC on board ROSAT 
and identified spectroscopically using the AUTOFIB fibre- 
optic system at the Anglo- Australian Telescope (A AT). Full 
details of the X-ray and optical observations are presented 
elsewhere (Georgantopoulos et al. 1994; Shanks et al. in 
preparation), and thus only brief details will be given here. 
The X-ray sample is based on 5 deep PSPC exposures of 
fields originally surveyed spectroscopically for QSOs as part 
of the Durham-UVX QSO survey (Boyle et al. 1990). The 
names, centres and exposure times for each PSPC field are 
given in Table 1. The position of the deepest field, GSGP4, 
is offset by 10 arcmin from the position of the SGP4 field 
centre in the Durham-UVX QSO sample. This was done to 
allow the GSGP4 field to overlap with a deep galaxy redshift 
survey (Broadhurst, Ellis & Shanks 1988) also made in this 
region. 

Over all 5 fields, 194 X-ray sources were detected at the 
5(7 limit in the 0.5-2 keV band. We use this band (rather 
than the full 0.1-2 keV band of the PSPC) to minimize any 
contribution from Galactic emission which can dominate be- 
low 0.5 keV. Due to the rapid increase in the size of the point- 
spread function (PSF) with off-axis angle, we limited our X- 
ray source detection to within 18 arcmin of the field centre in 
each PSPC image. This ensured that the PSF full width at 
half maximum intensity (FWHM) for any image is always 
less than two 15-arcsec pixels, and also avoided any sig- 
nificant obscuration by the telescope rib support structure 
which occurs at larger off-axis angles in the PSPC images. 
The total area surveyed at the 5(7 limit is therefore equal 
to 1.41 deg^. The flux limits corresponding to this detection 
limit not only vary from field-to-field (due to the different ex- 
posure times) but also vary within each field (due to the in- 
crease in the FWHM of the PSF with off-axis angle). To cal- 
culate the area covered by this survey as a function of X-ray 
flux, we therefore divided each field into 5 concentric annuli 
{6 < 10 arcmin, 10 < 61 < 12 arcmin, 12 < 6 < 14 arcmin. 



Table 1. ROSAT FSFC fields. 



Field 


RA (J2000) Dec 


Exposure Time 








(sees) 


SGP2 


00 52 04.8 


-29 05 24 


24494 


SGP3 


00 55 00.0 


-28 19 48 


21062 


GSGP4 


00 57 28.7 


-27 38 24 


48955 


QSFl 


03 42 09.6 


-44 54 36 


26155 


QSF3 


03 42 14.3 


-44 07 48 


27358 


14 < 61 < 16 


arcmin and 


16 < 61 < IS 


! arcmin) and deter- 



mined the flux limit in each annulus based on the field ex- 
posure time and mean theoretical FWHM of the PSF at this 
point. The resulting area coverage of the survey as a func- 
tion of X-ray flux limit is given by the first two columns of 
Table 2. 

2.2 Optical observations 

COSMOS and APM measurements of the UK Schmidt J 
plates in the South Galactic Pole region (J9771) and Field 
249 (J2762) were used to identify the optical counterparts to 
the X-ray sources in the survey. A transformation between 
the ROSAT X-iay and COSMOS/APM optical co-ordinate 
frames was achieved using the positions of the existing 10- 
12 Durham-UVX survey QSOs detected by ROSAT on each 
field (see Shanks et al. 1991). Optical counterparts to each 
X-ray source were then identified out to a radius of 45 arcsec 
from the transformed X-ray position. At the plate limit 
(B < 22.5 mag), the majority of X-ray sources had at least 
one optical counterpart within 20 arcsec. 

Low-resolution (12 A) optical spectra were obtained for 
the nearest optical counterpart to each X-ray source us- 
ing the AUTOFIB multi-object fibre-optic system at the 
AAT. In addition, we also obtained a few spectra for some 
of the fainter optical counterparts using the Low Disper- 
sion Survey Spectrograph (LDSS) also at the AAT. In 5 
nights we obtained optical identifications for the counter- 
parts to 145 of the 194 X-ray sources. Of the 49 sources 
with no optical identification, in 30 cases the spectra were 
too poor to permit a reliable identification. For the re- 
maining 19 objects, restrictions on the positioning of the 
fibres prevented the observation of the optical counterpart. 
The spectroscopic incompleteness was not independent of 
X-ray flux, ranging from < 10 per cent at 5(0.5 — 2 keV) > 
4 X lO'^^'ergs'^cm"^ to ~ 30 per cent at 5(0.5 - 2 keV) > 
4 X 10""'^® ergs~"'^cm~^. In the analysis below, we have cor- 
rected for this incompleteness by using an 'effective' survey 
area at each flux limit. For each flux limit in the survey, this 
'effective area' was simply obtained by multiplying the total 
survey area by the fraction of sources successfully identified 
to that flux limit. These effective survey areas are also listed 
in Table 2. 

Of the 145 sources with a reliable spectroscopic iden- 
tification, 107 (74 per cent) were classified as QSOs. These 
107 QSOs thus comprise the ROSAT QSO sample which will 
form the basis for the analysis presented below. An object 
was classified as a QSO if one or more broad emission lines 
(FWHM > 1000 km s~^) were present in the optical spec- 
trum. Based on the observed emission lines, we were able 
to determine an unambiguous redshift for 85 per cent of 
the QSOs. A further 12 galaxies exhibiting narrow emission 
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Table 2. Cumulative area coverage of ROSAT survey. 



Flux Limit 
5(0.5 - 2keV)ergs~lc 



Total area 
(deg2) 



Effective area 
(deg2) 



0.32> 


<10-14 


Q 


09 


0.06 


0.33> 


<10-14 





13 


0.09 


0.49> 


<10-14 





17 


0.12 


0.51> 


<10-14 





26 


0.18 


0.52> 


<10-14 





30 


0.21 


0.54> 


<10-14 





38 


0.27 


0.56X10-14 





42 


0.30 


0.60> 


<10-14 





48 


0.35 


0.63> 


<10-14 


Q 


56 


0.43 


0.64> 


<10-14 





60 


0.45 


0.66> 


<10-14 





69 


0.52 


0.68> 


<10-14 





73 


0.56 


0.75 > 


<10-14 





77 


0.61 


0.79> 


<10-14 





82 


0.64 


0.92> 


<10-14 


1 


02 


0.82 


0.10> 


<io-i3 


1 


07 


0.84 


0.11> 


<io-i3 


1 


12 


0.88 


0.12> 


<io-i3 


1 


18 


0.95 


0.14> 


<io-i3 


1 


24 


1.04 


0.15> 


<io-i3 


1 


30 


1.09 


0.17> 


<io-i3 


1 


36 


1.17 


0.18> 


<io-i3 


1 


41 


1.22 


0.30> 


<io-i3 


1 


41 


1.24 


0.40> 


<io-i3 


1 


41 


1.31 


0.65 > 


<io-i3 


1 


41 


1.41 



lines (FWHM < lOOOkms-l) were also identified in this 
survey. Without spectral coverage of the redshifted Ha/[Nll] 
region (and in some cases also the redshifted H/9/[OlIl] re- 
gion) it was difficult to determine the true nature of these 
galaxies (i.e. starbursts, Seyfert 2), although their strong 
stellar continua and weak emission lines are very similar to 
the properties of 'normal' late-type galaxies commonly seen 
in field galaxies redshift surveys at conducted at similar op- 
tical magnitude limits (B ~ 21 mag, see Broadhurst et al. 
1988). A more detailed discussion of their properties is pre- 
sented in Georgantopoulos et al. (1994), although the effect 
of including these galaxies in the ROSAT QSO sample is 
considered explicitly in section 3.2.2. 

As in Paper 1, we convert the 0.5-2 keV QSO fluxes to 
the Einstein 0.3-3.5 keV band using the relation: 

5(0.3 - 3.5 keV) = 1.85(0.5 - 2 keV) 

which is accurate to ±2 per cent for all X-ray spectral indices 
0.6 < ax < 1-5, fv oc . This enables us to combine our 

sample of 107 QSOs with the sample of 448 QSOs identified 
in the Emstem EMSS by Stocke et al. (1991). The EMSS 
QSO sample includes the 21 additional 'expected' QSOs in 
the EMSS, introduced to correct for the effects of incom- 
pleteness. In contrast to Paper 1, we choose to incorporate 
these QSOs in the analysis below in order to maintain con- 
sistency with the ROSAT QSO sample, which has also been 
corrected for incompleteness (albeit in a different fashion). 
We also note that the EMSS QSO sample also includes 31 
ambiguous and uncertain QSOs (tables 8 and 10 in Stocke 
et al. 1991). These objects exhibit predominantly narrow 
emission lines in their optical spectra, but were included in 
the QSO sample either because of suspected broad emission 
in the Balmer lines or because of high ionization implied 
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Figure 1. The number-redshift, n[z), diagram for the 107 QSOs 
identified in the ROSAT sample. The dashed and dotted lines 
represent the n[z) predictions for models H and K respectively 
(see text). 

from the [OlIl]/[Oll] line ratios. Some of these objects may 
thus be similar to the narrow emission line galaxies which 
we have also identified in our survey but chosen to exclude 
from our QSO sample. In the analysis below we will there- 
fore consider the EMSS QSO sample both with and without 
the inclusion of these ambiguous QSOs in order to ascertain 
the importance of these objects to the determination of the 
XLF and its evolution. 

In Fig. 1 we present the number-redshift, n{z), relation 
for the ROSAT sample. The n{z) relation exhibits a broad 
maximum over the redshift range 0.5 < z < 1.8. The me- 
dian redshift for the sample is z = 1.5, greater than the me- 
dian redshift for the EMSS, z = 0.2. The X-ray luminosity- 
redshift (Tx, z) diagram for the ROSAT and EMSS QSOs 
is shown in Fig. 2. The QSO X-ray luminosities have been 
calculated for go = 0-5 and Hq = 50kms-lMpc-l in the 
rest frame 0.3-3.5 keV energy band, assuming a power-law 
spectrum, fv (x with spectral index ax = 1- This is 

the mean spectral index for radio-quiet QSOs in the EMSS 
(Wilkes & Elvis 1987). Although the ROSAT QSOs ap- 
pear to exhibit a slightly softer mean spectrum, ax = 1.2 
(Stewart et al., in preparation), in order to maintain consis- 
tency with previous analyses (Paper 1; M91; Delia Ceca et 
al. 1992), we retain the assumption that ax = 1 throughout 
this paper. For the power-law luminosity evolution models 
considered below (ix oc (1 -|- z)^ , see Section 3.2.1), an in- 
crease in adopted value of ax simply increases the derived 
rate of evolution, k, by the same amount. 

In Paper 1 we also considered a two-component power- 
law model for the X-ray spectra of QSOs, which introduced 
significant softening (an increase in the value of ax) of the 
spectrum below 1.5 keV. With additional data and improved 
spectral fitting it is now apparent that, although the X-ray 
spectra of QSOs do indeed soften at low energies, this soft- 
ening occurs below 0.5 keV in the rest-frame of the QSOs 
(Stewart et al., in preparation), and is thus unimportant 
for the calculation of rest-frame luminosities in this analy- 
sis. From Fig. 2, we see that the addition of the ROSAT 
sample extends coverage of the QSO (ix) z) plane by more 
than an order of magnitude in luminosity. Although the 
ROSAT ssLiaple provides little improvement in the statistics 
of low-redshift (z < 0.2), low-luminosity (ix < 1043ergs-l) 



4 B.J. Boyle et al. 
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Figure 2. The X-ray luminosity-redshift diagram for the QSOs 
identified in the ROSAT sample (filled circles) and EMSS (open 
circles). 

QSOs, it contains 80 QSOs with z > 1, almost twice as many 
QSOs as the EMSS (46 QSOs) over the same redshift range. 

3 ANALYSIS 

3.1 1/Va analysis 

In order to obtain an initial estimate of the overall behaviour 
of the XLF and its dependence on redshift, we have derived 
binned estimates for the QSO XLF using the 1/Va statistic 
(Avni & Bahcall 1980). This method is fully described by 
M91 and Paper I. In Figs 3(a) and 3(b) we plot the derived 
0.3-3.5 keV QSO XLF for the combined ROSAT and EMSS 
(including ambiguous QSOs) samples for go = (Fig. 3a) 
and go = 0-5 (fig 3b). The corresponding QSO XLFs based 
on the ROSAT sample alone are shown in Figs 3(c) and 
3(d). In these figures, the XLFs are plotted in four redshift 
intervals of equal width in Iog(l + z) over the range Q < z < 
3. Error bars on these estimates of the XLF are obtained 
from Poissonian statistics. In this paper, we choose to quote 
the binned XLF per logarithmic luminosity interval rather 
than per linear luminosity interval adopted in Paper I and 
M91. The advantage of the approach adopted here is that 
luminosity evolution, in which the total number of QSOs 
is conserved with cosmic epoch, is represented by a purely 
horizontal shift of the XLF with redshift in this diagram. 
To avoid any misleading impression of the XLF coverage 
provided by these samples, all bins that contain three QSOs 
or less have been excluded from these figures. 

From Figs 3(a) and 3(b), we see that the XLF exhibits 
a characteristic two-power-law form, first identified by M91 
for the low-redshift (z < 0.2) XLF. In this case, however, 
the inclusion of the ROSAT sample has extended the XLF 
to lower luminosities at z > 0.2 and revealed that this two- 
power-law form is also apparent at these redshifts. This 
confirms the initial result reported in Paper I. The contri- 
bution of the ROSAT sample to the XLF is also borne out 
by Figs 3(c) and 3(d), where it can be seen that most of 



the QSOs identified in the ROSAT sample lie around the 
'break' luminosity or fainter for 0.4 < z < 2.9. From these 
figures it is also apparent that the redshift dependence of 
the XLF is governed predominantly by strong luminosity 
evolution for 2^ 1.8. This strong evolution is confirmed by 
the < VelVa > statistic (Avni & Bahcall 1980) calculated 
both for the combined data-set (< Ve/Va >= 0.67 ± 0.01) 
and for the /SO^AT sample alone (< Ve/Va >= 0.63 ±0.03). 
Both values are significantly different from the no-evolution 
value, < Ve/Va >= 0.5, at greater than the 4(7 level. The 
approximately constant luminosity shift between successive 
Alog(l -|- z) redshift bins at 2 < 1.8 suggests that the evo- 
lution could take a power-law form, ix c>c (1 -|- z)^ , at 
these redshifts. However, care should be taken not to over- 
interpret any binned representation of the XLF, as a sig- 
nificant amount of evolution can occur within each redshift 
interval plotted in such diagrams. Moreover, for luminos- 
ity evolution, even the higher value for < Ve/Va > (which 
primarily tests for density evolution) obtained for the com- 
bined sample compared to that derived for the /SO^AT sam- 
ple should not be taken to imply that the ROSAT sample 
exhibits less evolution than the EMSS. Indeed, the lower 
value simply reflects the fact that the ROSAT sample lies 
mostly on the flatter part of the XLF. 

3.2 Maximum likelihood analysis 

3.2.1 Models 

In order to derive a more accurate parametric representation 
of the QSO XLF and its evolution, we turn to the technique 
of maximum likelihood analysis employed in Paper I. This 
method, developed by Marshall et al. (1984) for use with the 
QSO optical luminosity function, yields best-fit values for 
parameters in a functional model used to describe the XLF 
and its evolution. As in Paper I, we use the two-power-law 
form for the z = XLF proposed by M91 (with modified 
normalization as described in Paper I): 

$x(ix) = $xix« < ix(^ = 0) 

^x(Xx)= ^,37.) X-J: ix>i^(^ = 0) 

X44 

where $x the normalization of the XLF and 71, 72 are 
the faint- and bright-end slopes respectively of the XLF. 
LX44 is 0.3-3.5 keV X-ray luminosity expressed in units of 
10'*'* erg s~^ . The evolution of the XLF is fully described by 
the redshift dependence of the 'break' luminosity, L^(z). As 
in Paper I, we investigate the two most commonly used func- 
tional forms to describe the evolution: power-law evolution 
in (l + z), 

ixW = ix(0) (1 + .)', 

and exponential evolution with fractional look-back time 
(r), 

L*^{z) = i^(0) exp(fcr). 

We also consider power-law evolution models in which the 
evolution 'switches off' at some maximum redshift Zmax- 

-^x(^) = ix(^max) Z > Zraax- 
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Figure 3. Binned 1/Va estimates for the QSO X-ray luminosity function in the 0.3— 3.5keV band based on: (a) EMSS and ROSAT 
samples for gg = O-O? (b) as (a) for gg = O-^? (c) ROSAT sample only for gg = 0-0' cl) (c) for 90 = O-^- 



Beyond Zmax the XLF is independent of redshift and 
the comoving space density of QSOs remains constant for 
any given X-ray luminosity. To complete the possible 
parametrizations of pure luminosity evolution, we also in- 
vestigate a 'polynomial' evolution model, 

similar to that adopted by Dunlop & Peacock (1990) in their 
study of the QSO radio luminosity function. 

In addition to pure luminosity evolution models, we 
have also investigated more complex evolutionary forms. We 
have tried to fit a model which incorporates both luminosity 
and density evolution. In this case the redshift dependence 
of the XLF is not only modelled by a (1 -|- z) power-law 
luminosity evolution (as above) but also includes sl (1 + z) 
power-law evolution in comoving number density: 

$xW = $x(0) (1 + .)*^=. 

We have also studied models in which the shape of the XLF 
is dependent on redshift, and have considered a power-law 



luminosity evolution model in which the faint-end slope also 
varies linearly with z: 

71(2) = 71(0) + ^22. 

Any potentially large dispersion in the X-ray spectral indices 
for QSOs could also affect the results derived in our analysis 
below. Giallongo & Vagnetti (1992) have demonstrated that 
any dispersion in the adopted spectral index for QSOs can 
also significantly reduce the derived rate of evolution. To 
estimate the size of this effect in the X-ray samples consid- 
ered below, we have included fits in which the dispersion in 
X-ray spectral index, a{ax), is assumed to be cr{ax) = 0.25 
and cr{ax) = 0.5. In order to calculate the best-fit model in 
these cases, we have used the expression derived by Francis 
(1993) to calculate the mean observed X-ray spectral index, 
aXoba ) ^ function of redshift introduced as a result of the 
dispersion in the spectral index: 

"X„b, = ax + o-^ {1 - 72) ln(l + z) . 

For an XLF steeper than 72 = 1 (as is observed here), any 
dispersion in the spectral index will therefore cause an ef- 
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fective 'hardening' in the mean observed spectral index with 
redshift, resulting in a corresponding decrease in the value of 
the evolution rate k. Although, this expression only strictly 
applies to a single-power-law LF, most QSOs in this anal- 
ysis lie on the steep portion of the XLF, and so the value 
of 72 can used for the slope without any significant loss of 
accuracy. Note that dispersion in spectral index does not 
affect the shape of the XLF and so the value of 72 used in 
this expression can be obtained from the value derived in 
the absence of any dispersion in the spectral index, with- 
out requiring any iteration of the above equation to obtain 
consistent values of «x„b, and 72 in the model fits. 

In the model fits to the XLF there can thus be four, 
five or six free parameters: 71, 72, L'^(Q), k, plus Zmax 
and/or k'j, depending on the evolutionary model chosen. 
The best-fit values for these parameters may be obtained 
from the maximum likelihood analysis. Statistical errors on 
these best-fit values (corresponding to 68 per cent confidence 
regions) were derived using the method described by Boyle, 
Shanks & Peterson (1988). However, the range in the values 
(particularly those for the z = XLF) from model to model 
may give a better indication of the full uncertainty in these 
parameters. The normalization of the XLF ($x) ^ 
free parameter in the fit and its value is derived from the 
normalization of the best-fit model to the total number of 
QSOs observed in the samples (see Paper 1). 

Although the maximum likelihood analysis gives best-fit 
values for the model parameters, it does not give a 'goodness 
of fit' for the best-fit model. To test the overall acceptabil- 
ity of the model fit to the data, we therefore used the two- 
dimensional (2D) KS statistic (Peacock 1985) employing the 
algorithm devised by Press et al. (1992). The KS probabil- 
ity (-Pks) is derived from comparison of the two-dimensional 
cumulative probability distributions in luminosity and red- 
shift for both the model and data. As in Paper 1, these dis- 
tributions were tested over the full luminosity and redshift 
range of the combined sample, 10'*^ < Lx < 10'*®ergs~^, 
< 2 < 3. 

3.2.2 Results 

The results of the maximum likelihood analysis and 2D KS 
tests applied to the combined EMSS and ROSAT samples 
are presented in Table 3 for a variety (a veritable A-Z) of 
the models discussed above. With the inclusion of the larger 
ROSAT sample in this paper, we now find that neither of 
the simplest evolutionary forms (exponential or power-law 
without a redshift cut-off) provides an acceptable fit (at the 
99 per cent confidence level) to the combined EMSS and 
ROSAT samples for either qo = or qo = 0.5 (models A- 
D). However, the KS probabilities of both the exponential 
and power-law evolution models are higher for the qo = 
case than for go = 0-5, in agreement with the trend found 
in Paper 1. Allowance for a dispersion in the X-ray spectral 
index for QSOs does not increase the acceptability of these 
simple models. Indeed, the inclusion of an increasingly large 
dispersion in the spectral index for the go = 0-5 power- 
law evolution model (models E and F) decreases the overall 
acceptability of the fit. 

Similarly, the inclusion of a cut-off in the evolution at 
high redshift does not result in an acceptable model fit to the 
data. For go = 0, the addition of a redshift cut-off (model 



G) only results in a very small increase in the overall KS 
probability of the model fit (from 0.5 per cent to 0.7 per 
cent). In contrast, for go = 0.5, the inclusion of a redshift 
cut-off in the power-law evolution model does significantly 
increase the acceptability of the model fit, from 0.02 per 
cent (model C) to 0.5 per cent (model H). Nevertheless, at 
this KS probability, this model can hardly be considered a 
good fit to the data. A better fit can be achieved if the X- 
ray spectral index is allowed to have a significant dispersion 
(models 1 and J). For a cr{ax) = 0.5, the KS probability for 
the go = 0.5, power-law evolution (k = 2.51) with redshift 
cut-off (zmax = 1.25) is 1.3 per cent, formally acceptable at 
the 1 per cent level. 

A similar level of acceptability, without the inclusion 
of a dispersion in the spectral index, is found for the poly- 
nomial evolution model in a go = universe (model K), 
although the corresponding model for go = 0.5 (model L) 
has a much lower KS probability (0.3 per cent). Inclusion of 
additional coefficients in the polynomial model (i.e. a third- 
order fit) yields no improvement in the overall acceptability 
of the model for go = 0.5. Indeed the formal best-fit value 
for the z^ coefficient in a cubic polynomial fit to the data in 
a go = 0.5 universe is ^3 = 0. It is intriguing that the val- 
ues for the evolution parameters k and k'j in the polynomial 
evolution model, k = 1.14(1.15) and k2 = -0.23(-0.26) 
for go = 0.0(0.5), are remarkably similar to those found by 
Dunlop & Peacock (1990), k = 1.18, k2 = -0.28 for the ra- 
dio evolution of flat-spectrum radio sources (predominantly 
QSOs) in a go = 0.5 universe. 

Given the poor fit of many of the luminosity evolution 
models applied to the data, we also considered the more 
complicated evolutionary models discussed above. However, 
as can be seen from Table 3, neither the luminosity and 
density evolution models (models M and N) nor the models 
incorporating a variable faint end slope (models O and P) 
provide improved fits to the combined data-set in either go = 
or go = 0.5 universes. 

Although we have accounted for X-ray flux-dependent 
incompleteness in the ROSAT survey (see Section 2.2), it 
is possible that some optical magnitude-dependent incom- 
pleteness has also been introduced by the decreasing suc- 
cess rate for optical identifications at faint magnitudes in 
the spectroscopic survey. Our correction for incompleteness 
(i.e. effective areas) assumes that these unidentified sources 
contain the same relative numbers of source types (QSOs, 
galaxies, etc.) as the identified sources. Since QSOs are 
generally the easiest class of object to identify in low signal- 
to-noise ratio spectra, this assumption may not be strictly 
correct. Currently 49 out of 194 objects remain unidentified, 
of which 30 are as a result of poor optical spectra, corre- 
sponding to spectroscopic incompleteness of 17 per cent. In 
order to establish whether any possible spectroscopic bias in 
this population significantly influences our results, we have 
therefore carried out the following two tests. 

First, we increased the X-ray flux limit in the ROSAT 
sample to 5(0.3 - 3.5 keV) > 2.0 x lO"^'* erg s"^cm"^, cor- 
responding to the point at which the spectroscopic incom- 
pleteness falls to below 10 per cent. At this limit, the 
ROSAT sample contains 105 X-ray sources (including 71 
QSOs), with 9 objects unidentified due to poor optical spec- 
tra and a further 10 sources not observed due to fibre posi- 
tioning restrictions. The results of performing the maximum 
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Table 3. QSO XLF models. 



Model 


go 


Evolution 




71 


72 


logL* (0)t 


k 


^max 




X + 


Pks 




A 


0.0 


(1 + z)^ 


0.00 


1.58 


3.35 


43.81 


2.63 






0.72 


4.8 X 


10-3 


B 


0.0 




0.00 


1.39 


3.30 


43.66 


5.10 






0.74 


3.5 X 


10-3 


C 


0.5 


(1 + zY 


0.00 


1.12 


3.31 


43.29 


2.44 






1.35 


1.6 X 


10-4 


D 


0.5 




0.00 


1.14 


3.33 


43.27 


4.02 






2.34 


1.7 X 


10-3 


E 


0.5 


(1 + zY 


0.25 


1.41 


3.29 


43.70 


2.20 






1.30 


2.6 X 


10-5 


F 


0.5 


(1 + zY 


0.50 


1.44 


3.29 


43.79 


1.54 






1.24 


4.3 X 


10-7 


G 


0.0 


(1 + zY + ^^max 


0.00 


1.53 


3.38 


43.70 


3.03 


1.89 




0.79 


6.7 X 


10-3 


H 


0.5 


(1 + zY + ^max 


0.00 


1.36 


3.37 


43.57 


2.90 


1.73 




1.45 


4.8 X 


10-3 


H' 


0.5 


(1 + zY + ^^max 


0.00 


1.37 


3.53 


43.54 


3.11 


1.83 




1.64 


2.7 X 


10-3 


H" 


0.5 


(1 + zY + ^max 


0.00 


1.36 


3.32 


43.59 


2.80 


1.75 




1.40 


3.2 X 


10-3 


I 


0.5 


(1 + zY + ^^max 


0.25 


1.45 


3.44 


43.61 


2.93 


1.45 




1.12 


5.4 X 


10-3 


J 


0.5 


(1 + zY + ^^max 


0.50 


1.46 


3.45 


43.65 


2.51 


1.25 




1.12 


1.3 X 


10-2 


K 


0.0 


kz + k2z'^ 


0.00 


1.50 


3.35 


43.71 


1.14 




-0.23 


0.84 


1.5 X 


10-2 


L 


0.5 


kz + k2z'^ 


0.00 


1.26 


3.32 


43.52 


1.15 




-0.26 


1.89 


2.6 X 


10-3 


M 


0.0 


(1 + zY + (l + zY^ 


0.00 


1.53 


3.35 


43.78 


2.76 




-0.26 


0.82 


4.5 X 


10-3 


N 


0.5 


(1 + zY + (1 + zY^ 


0.00 


1.55 


3.34 


43.81 


1.91 




1.03 


0.78 


1.2 X 


10-3 


O 


0.0 


(1 + zY + ^max + k2Z 


0.00 


1.46 


3.39 


43.65 


3.16 


1.92 


-0.18 


0.98 


5.8 X 


10-3 


P 


0.5 


(1 + zY + ^max + k2Z 


0.00 


1.39 


3.39 


43.63 


2.77 


1.60 


0.25 


1.16 


2.9 X 


10-3 


Q§ 


0.0 


(1 + zY 


0.00 


1.58 


3.35 


43.81 


2.66 






0.74 


2.2 X 


10-2 




0.5 


(i + zY 


0.00 


1.12 


3.31 


43.29 


3.17 






2.59 


1.4 X 


10-5 


s§ 


0.0 


(1 + zY + ^:max 


0.00 


1.53 


3.38 


43.70 


3.34 


1.79 




0.63 


3.1 X 


10-1 


T§ 


0.5 


(1 + zY + Zmax 


0.00 


1.36 


3.37 


43.57 


3.25 


1.60 




1.59 


1.2 X 


10-1 


UH 


0.5 


(1 + zY + Zmax 


0.00 


1.25 


3.40 


43.52 


3.17 


1.48 




1.61 


1.9 X 


10-3 




0.5 


(1 + zY + Zmax 


0.00 


1.59 


3.44 


43.63 


2.92 


1.75 




0.96 


1.2 X 


10-2 


w+ 


0.0 


(1 + zY + Zmax 


0.00 


1.53 


3.18 


43.81 


2.14 


1.25 




0.76 


6.3 X 


10-2 


x+ 


0.5 


(1 + zY + Zmax 


0.00 


1.02 


3.02 


43.41 


2.49 


0.95 




3.39 


1.0 X 


10-2 




0.0 


(1 + zY + Zmax 


0.00 


1.60 


3.51 


43.88 


2.93 


1.80 




0.95 


2.0 X 


10-2 


z* 


0.5 


(1 + zY + Zmax 


0.00 


1.57 


3.81 


43.79 


3.13 


1.58 




1.27 


5.9 X 


10-2 


Errors 








±0.15 


±0.1 


±0.2 


±0.1 


±0.1 


±0.1 









t In units of ergs . 

t In units of 10-'^ Mpc-3(10'14 ergs- 1)- 1 . 

' ROSAT sample limited to 5(0.3 - 3.5keV) > 2 X lO-l'* ergs-lcm-2. 
" Excluding all unidentified sources in the ROSAT sample. 
%ROSAT sample only. Parameters for z = XLF held fixed in fit. 
^ Excluding EMSS 'ambiguous' QSOs. 

Including ROSAT narrow emission line galaxies. 
~^ No correction applied to ROSAT 0.5— 2keV fluxes. 

* EMSS sample limited to 5(0.3 - 3.5keV) > 2.8 X 10- 13 ergs- lcm-2 . 



likelihood analysis for the power-law evolution model with 
redshift cut-off on this restricted /SO^AT sample, combined 
with the EMSS, are listed in Table 3 (model H'). Compared 
with the corresponding model for the full ROSAT sample 
(model H) we can see that the use of the restricted sam- 
ple has made no significant difference (at > 2a level) to the 
overall model parameters, or to the acceptability of the fit. 

Secondly, we also performed the maximum likelihood 
analysis for the same evolution model on the full ROSAT 
sample, assuming that none of the unidentified sources was 
a QSO (model H"). Although this assumption is almost 
certainly unrealistic (since it also includes sources that were 
simply not observed), it does provide an upper limit in the 
more likely case that the majority of the unidentified sources 
are objects with weak emission line or absorption spectra 
(e.g. galaxies, clusters of galaxies, BL Lacs) and not QSOs. 
However, as can be seen from the values of best-fit param- 
eters for model H" listed in Table 3, even this pessimistic 
assumption makes no significant difference to the overall re- 
sult. From these two tests we therefore conclude that any 
optical magnitude effects introduced by the spectroscopic 



incompleteness do not significantly influence our results. 
4 DISCUSSION 

4.1 ROSAT versus EMSS evolution 

It is clear from the previous section that few of the models 
investigated in this paper provide an acceptable fit to the 
combined EMSS and /SO^AT samples. Although pure lumi- 
nosity evolution models can be found that will formally fit 
the data at the 1 per cent level for both qo = (polynomial 
evolution, model K) and go = 0-5 (power-law evolution with 
large dispersion in X-ray spectral index, model J), neither 
of the models is acceptable at the 5 per cent level. Indeed, 
more complex evolutionary forms (e.g. luminosity and den- 
sity evolution, models M and N) provide no better fits to 
the combined data-set. 

Nevertheless, simple evolution models do indeed pro- 
vide good fits (-Pks > 10 psr cent) to both EMSS and 
ROSAT data-sets when the maximum likelihood analysis is 
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performed separately on each sample. In Paper 1, we demon- 
strated that the EMSS QSOs were consistent with power-law 
evolution, L*y^{z) oc [l + zf ''^' for go = 0.5, with the ROSAT 
sample (which then contained only 42 QSOs) exhibiting a 
much higher rate of evolution, L^(z) oc (1 -|- z)^'^ . The rate 
of evolution was derived for the ROSAT sample by adopt- 
ing the best-fit values for the parameters of the z = XLF 
obtained for the combined EMSS and /SO^AT samples (71, 
72 and i^(0)) and deriving only the value of k from the 
maximum likelihood analysis. If we adopt a similar proce- 
dure for the larger ROSAT data-set, we again find that we 
derive much higher rates for the evolution of the XLF, with 
L'^(z) oc (1 + 2)2^^^ and L'^(z) oc (1 + zf '^'^ for go = 
(model Q) and go = 0.5 (model R) respectively. 

However, neither model is a good fit to the ROSAT 
data: the go = model is rejected at the 95 per cent con- 
fidence level, while the go = 0.5 model is rejected at the 99 
per cent confidence level. A much better fit to the ROSAT 
sample is obtained for models with a cut-off in the luminos- 
ity evolution at Zmax (models S and T). From the maximum 
likelihood analysis, we obtain best-fit parameter values of 
2max = 1.6(1.7) and k = 3.25(3.34) for go = 0.5(0.0). Both 
fits have a much higher KS probability (> 10 per cent) than 
any of the other models and confirm that pure luminosity 
evolution with a redshift cut-off is an acceptable fit to the 
ROSAT sample. 

The value for the evolution parameter (k = 3.34 ± 0.1) 
derived from the /SO^AT sample alone (go = 0.0) is signifi- 
cantly higher than that derived from the EMSS QSO sample 
(k = 2.56), or even from the combined EMSS/ ROSAT ssLia- 
ple (k = 3.03 ± 0.1). The high rate of evolution implied 
by the ROSAT sample is also apparent from the n{z) dia- 
gram plotted in Fig. 1, where the observed median redshift 
for the sample is higher than than predicted by the power- 
law evolution fit (model H, dashed line) to the combined 
data-set. However, a significantly better fit to the observed 
n{z) is achieved with the go = polynomial evolution model 
(model K, dotted line). 

In Fig. 1 we have normalized both the predicted 
number-redshift relations (models H and K) to the num- 
ber of QSOs observed in the ROSAT sample. However, it 
is a general feature of many of the models investigated here 
that they significantly under-predict the number of QSOs 
in the ROSAT sample. For example, the integral log N - 
log S, N{> S), relation predicted by model H (dashed line 
in Fig. 4) falls significantly below the N{> S) relation for 
the ROSAT sample (filled circles), while providing a good 
fit to the N{> S) relation for the EMSS (triangles). Based 
on model H, we would predict only 75 QSOs in the ROSAT 
sample, compared to the 107 observed. From Poisson statis- 
tics, this is a 3a discrepancy and is consistent with the low 
KS probability obtained for model H. However, low count 
predictions for the ROSAT sample are not a universal fea- 
ture of the models considered above. The polynomial evo- 
lution model K (represented by the dotted line in Fig. 4) 
provides a much better fit to the ROSAT counts, predict- 
ing 119 QSOs in the ROSAT sample, consistent with the 
number observed to within ~ la. 

From Fig. 1 it is apparent that much of the discrepancy 
between the samples for power-law evolution models occurs 
in the number of QSOs predicted/observed at low redshift. 
It is therefore possible that the differences in the classifi- 



cation of low-redshift (z < 0.5) objects with narrow emis- 
sion lines between the two surveys (see Section 2.2) could 
explain some of this discrepancy. To estimate the size of 
this effect, we performed the maximum likelihood analysis 
for the power-law evolution model with a redshift cut-off 
for the combined EMSS and /SO^AT samples, both exclud- 
ing the ambiguous EMSS sources (model U) and including 
the 12 emission-line galaxies identified in the /SO^AT sam- 
ple (model V) by Georgantopoulos et al. (1994). As can 
be seen from the reported fits in Table 3, the exclusion of 
the AGN classified as 'ambiguous' in the EMSS makes little 
difference to the overall acceptability of the fit (cf. model 
H). In contrast, the inclusion of the emission-line galaxies 
in the ROSAT sample does increase the KS probability of 
the fit from 0.0048 to 0.012. Despite the marginal accep- 
tance of the model under these circumstances, we consider 
it unlikely that the emission-line galaxies identified in the 
ROSAT sample are related to AGN. As discussed in Sec- 
tion 2.2, all members of this class have spectra characteris- 
tic of late-type galaxies (see Georgantopoulos et al. 1994), 
with strong stellar continua and narrow (< 500 km s"^) [On] 
emission lines whose equivalent widths (~ 20 A) are consis- 
tent with 'normal' emission-line galaxies identified in field 
galaxy surveys at similar optical magnitude limits (see e.g. 
Broadhurst et al. 1988). Clearly, better optical spectra will 
be required to determine the true nature of this population. 

The relative success of the polynomial luminosity evo- 
lution model over the power-law evolution model (with a 
redshift cut-off) suggests that the latter model is no longer 
an adequate representation of QSO evolution in the X-ray 
regime. Nevertheless, the power-law evolution model does 
provide a good fit to both the EMSS and ROSAT samples 
when they are analysed separately (see above). It is there- 
fore possible that the origin for the discrepancy lies not in 
the inadequacy of the model but in a systematic effect be- 
tween the EMSS and /SO^AT samples (possibly introduced 
by the different energy ranges for the two samples) which 
has not been accounted for in the analysis above. The most 
straightforward origin for any such effect lies in an error in 
the flux scales for one or both of the samples used in this 
analysis. However, a power-law luminosity evolution model 
can only be made compatible with both the EMSS and 
ROSAT samples if no correction is applied to the ROSAT 
0.5-2 keV fluxes in order to convert to the EMSS 0.3-3.5 keV 
band, although such models (W and X in Table 3) are only 
acceptable at the ^ 1 per cent level. Since even a radical al- 
teration in the assumed X-ray spectral index for QSOs can- 
not significantly alter the bandpass correction (see Boyle et 
al. 1993), such a result would require an ~ 80 per cent rela- 
tive shift between the Einstein and ROSAT Hnx scales (with 
the current ROSAT scale being too bright). While residual 
errors at less than the 20 per cent level may still occur in 
the flux scales, it is unlikely that an error of this magnitude 
is still present. Different spectral indices in the rest-frame 
spectra of the EMSS and ROSAT QSOs (at ~ 1.5 keV and 
~ 2.5 keV respectively) could also account for some of the 
difference in the derived evolution rates. Although this ex- 
planation was explored in Paper 1, it appears unlikely that 
the rest-frame X-ray spectra of QSOs harden sufficiently 
(Aax ~ — 0.5) at > 2 keV (Stewart et al. in preparation) 
for this to be a viable explanation. Alternatively, because 
the ROSAT QSOs are mostly observed at higher rest-frame 
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Figure 4. The log N- log S, N(> S), relations for the ROSAT 
sample (filled circles) and EMSS (triangles). The A^(> S) pre- 
dictions based on models H and K (see text) are denoted by the 
dashed and dotted lines respectively. 

energies than the EMSS QSOs, a larger dispersion in ax 
at harder energies would decrease the power-law evolution 
rate derived from the ROSAT sample, forcing it into better 
agreement with that obtained for the EMSS. The increase 
in the cr{ax) would, however, need to be dramatic. Based 
on the results presented in Table 3, the required decrease in 
the power-law evolution parameter to achieve consistency 
between the EMSS and ROSAT samples (Ak ~ 0.4) could 
only be obtained by effectively increasing the dispersion in 
the spectral index from cr{ax) = at 1.5 keV to cr{ax) = 0.5 
at 2.5 keV. 

It has also recently been claimed by Franceschini et al. 
(1994) that significant incompleteness exists in the EMSS 
at flux levels 5(0.3 - 3.5 keV) < 2.8 x 10"^^ erg s'^cm"^. 
By excluding QSOs in the EMSS with 5(0.3 - 3.5 keV) < 
2.8 X 10~-^^ergs~-^cm~^, Franceschini et al. (1994) were 
able to show that the remaining QSOs were consistent with 
a much higher rate of evolution Lx oc (1 -|- z)^'^ , for a spec- 
tral index ax = 1-2. This is equivalent to Lx oc (1 -|- z)^'^ 
for ax = 1.0 adopted here. We have carried out a max- 
imum likelihood analysis on the combined ROSAT and 
EMSS data-sets, with the EMSS restricted to flux levels 
5(0.3 - 3.5 keV) > 2.8 x 10"^^ erg s'^cm"^ . At this flux 
limit the EMSS contains 289 QSOs. The best-fit parameters 
for a power-law evolution model with a redshift cut-off ap- 
plied to this data-set are listed in Table 3 (models Y and Z) 
for both go = and go = 0-5. We find that the KS probabil- 
ities of the fits are significantly improved over the equivalent 
models (G and H) fitted to the entire EMSS, with the KS 
probability for the go = 0.5 model over ~ 5 per cent. The 
evolution rate derived for the go = 0.5 model, k = 3.13±0.1, 
is consistent with the value derived by Franceschini et al. 
(1994), although precise comparison is difficult as the model 
fitted here includes a redshift cut-off at Zmax- 

Finally, Hasinger et al. (1993) find some evidence for 
hardening of source spectra in their Lockman Hole survey at 
fluxes below 1.0 x 10"^'* erg s~^cm~^. Any hardening of the 
QSO spectra (and consequent decrease in the derived value 
of k) at such faint flux levels in the ROSAT sample could 



therefore help to explain some of the discrepancy between 
the EMSS and ROSAT evolution rates. However, we find 
no evidence for any significant increase in the QSO hardness 
ratios over the 0.5-2 keV band at similar flux levels in our 
survey (Almaini, private communication). It is therefore 
probable that sources other than QSOs are responsible for 
the hardening seen in the Hasinger et al. (1993) sample, 
which is currently largely unidentified. 

4.2 The 1 2 keV X-ray background 

We can also compute the expected contribution of QSOs 
to the X-ray background (XRB) based on the observed ex- 
tragalactic 1-2 keV background intensity, /xRB(1^2keV) = 
1.25 X 10~® erg cm~^s~"'^sr~"'^ , obtained by Hasinger et al. 
(1993). Despite the relatively poor fit of the XLF to the 
data, we may still obtain some useful limits based on the 
best of these models, since many of the parameters derived 
in this paper (particularly the evolution at high redshift) are 
more tightly constrained than in Paper 1. We compute the 
QSO contribution (/q) to the 1-2 keV XRB by integrating 
over the z = XLF as follows: 



ix(0)„ 



0.28 c 



ix$x(ix,0)(l-F2) 



-dzdL 



4T-ffo J J + 2qoz){l + 2)2 + «x 

ix(0)„i„ ^mi„ 

where the integration proceeds over the un-evolved luminos- 
ity range 10^^ < Lx{0) < lO'^^erg and over the redshift 
range < z < A. The factor 0.28 converts the background 
derived from the 0.3-3.5 keV band XLF to the 1-2 keV band, 
assuming a spectral index of ax = 1. 

Table 4 lists the values of /q obtained for five of the 
best-fit evolutionary models in this paper (models G, H, 
K, S and T). Beside each estimate of /q in this table, we 
also give the corresponding fractional contribution of QSOs 
to the 1-2 keV XRB. The QSO contribution ranges from 



0.42 X 10" 



0.66 X lO'^'keVcm'^s" 



or 34-53 per 

cent of the 1-2 keV XRB. The best-fit polynomial evolution 
model to the combined ROSAT and EMSS data-set (model 
K, go = 0) predicts a contribution of 50 per cent to the 1- 
2 keV XRB. The range is lower than that obtained in Paper 
1, partly because it has become clearer that a redshift cut-off 
is required in the evolution of QSOs at high redshifts, but 
principally because the faint-end slopes of the XLF derived 
in this paper are flatter than those obtained in Paper 1 (e.g. 
71 = 1.53 for go = in models G and S, cf. 71 = 1.71 
for the corresponding models in Paper 1). Nevertheless, the 
dominant uncertainty in the QSO contribution to the XRB 
is still the faint-end slope of the XLF. For example, a steep- 
ening of the XLF faint-end slope in model T by 0.15 (the 
1(7 error) to 71 = 1.51, whilst retaining the same values for 
the other parameters, would result in a prediction of 62 per 
cent for the QSO contribution to the XRB. Conversely, a 
flattening of the slope by the same amount to 71 = 1.21 
would give only 33 per cent. In contrast, an increase in the 
maximum redshift for the integration of the XLF from z = 3 
to 2 = 5 in model T only results in an increase for the QSO 
contribution to the XRB from 42 per cent to 46 per cent. 

A recent observation of the XRB with the ASCA satel- 
lite by Gendreau et al. (1994) yields a 1-2 keV back- 
ground, /xRB(l-2keV) ~ 1.21 X 10"® ergcm"^s"^sr"^ 



10 B.J. Boyle et al. 



Table 4. The QSO contribution to the 1-2 keV X-ray back- 
ground. 



Model 


(keV cm' 


-2s-lsr-l) 


/ ^XRB 


G 


0.64 


X 10-8 


0.52 


H 


0.42 


X 10-8 


0.34 


K 


0.63 


X 10-8 


0.50 


S 


0.66 


X 10-8 


0.53 


T 


0.56 


X 10-8 


0.44 



(based on a normalization of /xRB(lkeV) = 1.39 x 
10-8 erg cm-^s-^sr-^, ^ spectral index ax = —0.35, 

see table 1 in Gendreau et al. 1994), only 4 per cent less 
than that derived by Hasinger et al. (1993). This is de- 
spite the significantly flatter slope (j/-*^^®) measured for 
the XRB in this energy range by Gendreau et al. (1994). 
Based on these observations, the QSO contribution to the 
1-2 keV background would be correspondingly increased by 

4 per cent for all models listed in Table 4. 

5 CONCLUSIONS 

We have investigated the evolution of the XLF based on a 
new sample of 107 ROSATli.-Ta.y QSOs identified at faint X- 
ray flux levels, 5(0.5 - 2 keV) > 4 x 10-^® ergs-^cm-^. For 
go = 0-5, the evolution of the sample is consistent at the 95 
per cent confidence level with power-law luminosity evolu- 
tion, L'^{z) oc (l-l-z)^ 25 at 2 < 1.6 and a constant comoving 
space density at higher redshifts. This is significantly higher 
than the rate of evolution derived previously for the EMSS. 
Indeed, with the exception of a polynomial evolution model 
{L*y^{z) oc io(l 14^-0 23^=)^ ^ ^ powerTaw evo- 

lution model (ix(-^) (1 -|- z)^'^, Zmax = 1.25, go = 0-5), 
which also requires a large dispersion in X-ray spectral in- 
dex, (T(ax) = 0-5, most pure luminosity evolution models do 
not provide a good fit to the combined ROSAT and EMSS 
QSO samples. More complex evolutionary forms (luminosity 
and density evolution, a redshift-dependent faint-end slope 
of the XLF) also fail to provide an adequate overall fit to 
the combined data-set. It is possible that systematic effects 
introduced by the different energy bands sampled by the 
EMSS and ROSAT surveys (e.g. an increase in the X-ray 
spectral index dispersion with energy) and/or possible in- 
completeness in the EMSS could account (at least in part) 
for the poor fit of many models to the combined EMSS and 
ROSAT data-set. Based on extrapolation of the best-fit 
models for the XLF and its evolution, the most likely val- 
ues for the QSO contribution to the 1-2 keV XRB lie in the 
range 34-53 per cent. 
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